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HIGHLIGHTS 


►  Oxygen  dissolution  process  from  gas  phase  into  the  ionomer  was  included  in  the  model. 

►  Dissolution  kinetics  were  assessed  using  thin  ionomer  films  on  a  Pt  electrode. 

►  Agglomerate  sizes  were  determined  based  on  diffusion-limited  currents. 

►  The  estimated  agglomerate  size  is  close  to  that  observed  in  the  catalyst  layer. 

►  Slow  dissolution  is  the  main  hindrance  to  oxygen  transport. 


ARTICLE  INFO 


ABSTRACT 


Article  history: 

Received  7  May  2012 
Received  in  revised  form 
3  August  2012 
Accepted  22  August  2012 
Available  online  10  September  2012 


Keywords: 

Polymer  electrolyte  fuel  cell 
Oxygen  transport 
Catalyst  layer 
Agglomerate  size 
Dissolution  rate 


A  model  for  oxygen  transport  in  the  catalyst  layer  of  polymer  electrolyte  fuel  cells  was  developed.  The 
model  includes  oxygen  transport  resistance  for  dissolution  from  the  gas  phase  into  the  ionomer,  which  is 
typically  assumed  to  be  negligible  in  conventional  models.  The  dissolution  kinetics  were  experimentally 
assessed  by  measuring  the  diffusion-limited  current  density,  id,  in  a  planar  Pt  electrode  covered  with  thin 
ionomer  films  as  a  function  of  film  thickness.  The  extrapolation  of  i^1  to  zero  thickness  reflects  the 
dissolution  resistance.  The  effect  of  dissolution  resistance  is  added  to  a  conventional  agglomerate  model 
of  the  catalyst  layer  as  an  additional  cause  of  transport  loss.  The  agglomerate  size  was  determined  so  that 
the  model  predicted  the  same  diffusion-limited  current  density  as  the  experimental  results.  The 
agglomerate  size  was  close  to  the  particle  size  observed  in  scanning  electron  micrographs  of  the  cross- 
section  of  the  catalyst  layer  prepared  by  focused  ion  beam  milling.  In  contrast,  the  conventional 
agglomerate  model  must  adopt  an  agglomerate  size  that  is  larger  by  one  order-of-magnitude  than  that 
experimentally  observed.  As  a  result,  it  is  reasonable  to  attribute  the  hindrance  of  oxygen  transport  in  the 
catalyst  layer  to  the  slow  oxygen  dissolution  at  the  gas-ionomer  interface,  rather  than  to  the  presence  of 
large  agglomerates. 

©  2012  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

The  output  voltage  of  a  polymer  electrolyte  fuel  cell  (PEFC)  is 
strongly  affected  by  the  mass-transport  limitation  of  oxygen  at  high 
current  densities,  due  to  a  decrease  in  the  oxygen  concentration  at 
the  catalyst  surface.  The  oxygen  concentration  profile  has  been 
estimated  using  a  multicomponent  gas  diffusion  model  [1]  and 
Knudsen  diffusion  model  [2],  which  reflect  the  transport  of  species 
in  the  pores  of  the  electrodes.  Agglomerate  models  [3,4],  which  will 
be  reviewed  and  modified  in  this  paper,  have  been  employed  to 
calculate  the  oxygen  concentration  in  the  ionomer  that  covers  the 
catalyst.  A  clear  understanding  of  where  and  how  the  concentration 
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decreases  and  identification  of  the  rate-determining  process  in  the 
electrodes  will  facilitate  the  material  development  and  component 
optimization  of  PEFCs. 

Agglomerate  models  were  introduced  to  predict  current- 
voltage  characteristics  [3,4]  on  the  basis  of  scanning  electron 
microscopy  (SEM),  which  indicated  the  agglomeration  of  Pt/C 
particles  in  the  catalyst  layer.  The  agglomerate  model  ascribes  the 
increase  in  the  Tafel  slope  with  increasing  current  density  to 
a  decrease  in  the  oxygen  concentration  near  the  center  of  the 
agglomerate.  Mathematically,  the  agglomerate  is  simplified  into 
a  sphere  (or  a  cylinder  without  pores)  of  a  homogeneous  mixture 
consisting  of  a  carbon-supported  catalyst  and  an  ionomer.  The 
agglomerate  may  be  surrounded  by  a  film  of  the  ionomer.  When  the 
current  density  increases,  the  reaction  preferably  occurs  near  the 
surface  of  the  agglomerate  and  most  of  the  oxygen  is  consumed 
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there.  The  model  was  validated  by  a  comprehensive  comparison 
between  the  model  predictions  and  experimental  results  [5].  An 
agglomerate  radius  of  between  0.1  and  1  pm  was  typically  assigned 
in  the  previous  studies  [4-27]. 

If  this  model  is  sufficient  to  depict  the  catalyst  layer,  then  the 
same  size  agglomerates  should  be  observed  in  actual  catalyst 
layers,  and  consequently,  the  agglomerates  can  be  modified  to 
improve  the  catalyst  layer  according  to  the  model  predictions.  To 
our  knowledge,  no  such  study  has  been  conducted  and  the 
agglomerate  radius  has  been  treated  simply  as  a  parameter  that  fits 
the  model  predictions  to  the  experimental  results. 

This  work  proposes  a  model  for  the  oxygen  transport  in  the 
cathode  catalyst  layer.  Cross-sections  of  catalyst  layers  with  various 
ionomer  loadings  were  observed  to  estimate  the  agglomerate  size. 
The  oxygen  transport  resistance  of  the  ionomer  was  measured  for 
thin  films  on  an  electrochemical  cell  and  was  separated  into 
internal  and  external  resistances.  The  latter  is  incorporated  into  the 
conventional  agglomerate  model  as  an  additional  cause  of  trans¬ 
port  loss.  The  proposed  model  was  assessed  according  to  the 
current— voltage  characteristics  and  the  corresponding  agglom¬ 
erate  size. 

2.  Experimental 

2.1.  Fabrication  and  characterization  of  membrane-electrode 
assemblies 

Membrane-electrode  assemblies  (MEAs)  with  13  cm2  active 
areas  were  fabricated  by  hot  pressing  catalyst  layer  decals  onto 
25  pm  thick  Nation  membranes  (NRE-211,  DuPont).  The  anode 
(cathode)  decals  consisted  of  10  wt%  Pt/C  (TEC10E10A,  Tanaka 
Kikinzoku  Kogyo)  and  a  Nation  dispersion  solution  (DE2020  CS 
type,  DuPont),  and  30  wt%  Pt/C  (TEC10E30E,  Tanaka  Kikinzoku 
Kogyo)  and  DE2020  CS,  respectively.  Three  cathode  decals  with 
different  ionomer/carbon  weight  ratios  (I/C)  were  applied.  Plat¬ 
inum  loading,  I/C,  thickness,  and  the  porosity  of  the  resulting  anode 
and  cathode  catalyst  layers  are  listed  in  Table  1.  Developmental  gas 
diffusion  layers  were  inserted  when  the  MEA  was  assembled  into 
a  cell  fixture. 

The  cathode  decals  used  for  the  MEAs  in  Table  1  were  cut  out 
and  hot-pressed  onto  Nation  membrane  for  morphological  obser¬ 
vations.  Cross-sections  perpendicular  to  the  membrane  surface 
were  prepared  by  focused  ion  beam  (FIB;  FB200A,  Hitachi  High- 
Technologies)  milling.  Cross-sections  were  observed  using  SEM 
(S-5500,  Hitachi  High-Technologies). 

2.2.  Measurement  of  the  oxygen  transport  properties  of  thin 
ionomer  films 
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Fig- 1.  Schematic  illustration  of  the  cell  used  to  measure  oxygen  transport  properties  of 
an  ionomer  film,  (a)  Top  view  of  the  cell,  (b)  enlarged  view  of  the  electrodes,  (c)  side 
view  around  the  working  electrode  B,  and  (d)  experimental  setup. 


The  oxygen  transport  characteristics  were  measured  using  thin 
ionomer  films  prepared  on  an  electrochemical  cell  (011598,  BAS),  as 
shown  in  Fig.  1(a)  and  (b),  in  which  Pt  film  electrodes  were  printed 
on  a  quartz  plate.  Three  out  of  eight  electrodes  were  platinized  and 
used  for  the  measurement.  Electrodes  A  (1  mm  x  2.5  mm),  B 
(10  pm  x  2.5  mm),  and  C  (10  pm  x  2.5  mm)  served  as  the  counter, 


working,  and  pseudo-reference  electrodes,  respectively.  The  gap 
between  the  working  and  pseudo-reference  electrodes  was  30  pm. 
With  0.1%  oxygen,  the  quasi-steady-state  diffusion-controlled 
current  density  of  the  working  electrode  would  be  approximately 
1  A  cm-2  between  1  and  1000  s  using  the  value  of  gas  diffusivity 
[28].  Nation  films  were  spin-coated  on  the  electrodes,  as  shown  in 


Table  1 

Parameters  of  MEAs. 


MEA 

Anode 

Cathode 

Pt  loading  (mg  cm-2) 

I/C 

Thickness 

(pm) 

Porosity 

(%) 

Pt  loading  (mg  cm-2) 

I/C 

Thickness 

(pm) 

Porosity 

(%) 

Ionomer  volume  fraction 

(%) 

ECS  A  (cmpt  cm  2) 

a 

0.025 

1.0 

8 

71 

0.094 

0.5 

5 

66 

11 

42 

b 

0.025 

1.0 

8 

71 

0.094 

0.75 

5 

60 

17 

40 

c 

0.023 

1.0 

8 

74 

0.109 

1.0 

6 

56 

22 

49 
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Fig.  1(c),  at  room  temperature  using  Nation  dispersion  solution 
(DE2020,  DuPont)  diluted  with  1 -propanol.  Unwanted  film  in  the 
peripheral  area  was  rubbed  off  with  a  polyethylene  foam  swab 
dampened  with  ethanol.  The  film  was  further  dried  in  air  at  130  °C 
for  1  h.  Fig.  1(d)  shows  the  experimental  setup  used  for  the 
measurement.  The  cell  was  enclosed  in  a  glass  container  with 
a  nominal  capacity  of  200  mL  and  placed  in  a  temperature- 
controlled  oven.  Oxygen  was  diluted  with  nitrogen  to  1%  (for 
films  thicker  than  400  nm)  or  0.1%  (for  films  thinner  than  400  nm) 
and  supplied  to  the  cell  at  100  mL  min-1.  The  relative  humidity  of 
the  gas  was  controlled  by  mixing  fully  humidified  gas  and  dry  gas. 
Cyclic  voltammetry  (CV)  was  performed  at  5  mV  s-1  between  the 
potential  just  above  hydrogen  evolution  and  100  mV  vs.  the 
pseudo-reference  electrode  using  an  electrochemical  measurement 
system  (HZ-5000/HAG1512m,  Hokuto  Denko).  The  thickness  of  the 
films  was  measured  using  a  stylus  profiler  (Dektak  3  Surface 
Profiler,  Veeco)  after  the  electrochemical  measurement. 

For  comparison,  the  solubility  and  diffusion  coefficient  of 
oxygen  were  measured  for  a  96  pm  thick  cast  Nation  film  using 
a  microelectrode  cell.  Fig.  2(a)  shows  a  schematic  of  the  experi¬ 
mental  setup.  In  this  cell,  the  Nation  film  is  sandwiched  between 
a  50  pm  diameter  Pt  microelectrode  and  a  polytetrafluoroethylene 
(PTFE)  filter  through  which  oxygen  is  supplied.  The  edge  of  the  film 
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is  attached  to  NRE-211  Nafion  membrane  with  counter  (platinized 
Pt)  and  reference  electrodes  (RHE)  attached,  as  shown  in  Fig.  2(b). 
The  solubility  and  diffusion  coefficient  were  determined  from 
chronoamperometric  measurements  with  a  potential  step  applied 
from  1.1  to  0.4  V  vs.  RHE  [29].  The  measurement  was  conducted  at 
80  °C  under  three  different  relative  humidity  (RH)  conditions. 

2.3.  Measurement  of  the  MEA  performance 

Electrochemical  measurements  of  the  MEAs  were  conducted 
using  a  charge/discharge  unit  (5  V,  30  A,  Hokuto  Denko)  with  a  fuel 
cell  testing  system  (Chino).  The  reactant  gas  (hydrogen  or  oxygen) 
and  nitrogen  were  mixed  and  humidified  by  bubbling  through 
water  at  controlled  temperatures.  The  flow  rates  of  the  anode  and 
cathode  dry  gases  were  always  set  at  500  and  2000  standard  cubic 
centimeter  per  minute  (SCCM),  respectively,  to  ensure  reproduc¬ 
ibility  of  the  RH  and  a  sufficiently  high  stoichiometric  ratio.  A  high 
stoichiometric  ratio  helps  to  justify  the  omission  of  the  in-plane 
variation  in  the  mathematical  model  described  later.  The  partial 
pressures  of  the  gases  at  the  inlet  can  be  calculated  from  the  flow 
rates  of  the  gases,  the  water  vapor  pressure,  and  the  total  pressure. 

The  current— voltage  characteristics  of  the  MEA  and  the  elec¬ 
trochemical  surface  area  (ECS A)  of  the  cathode  were  measured  at 
80  °C  and  140  kPa.  Firstly,  hydrogen-nitrogen  and  oxygen- 
nitrogen  mixtures  with  70%  RH  were  supplied  to  the  anode  and 
cathode,  respectively.  The  RH  was  selected  to  ensure  that  the  cell 
was  free  from  flooding  over  the  entire  range  of  the  measured 
current  density.  The  partial  pressure  of  hydrogen  was  set  at  93  kPa, 
and  that  of  oxygen  was  set  at  1.9,  9.3,  or  19.5  kPa.  The  current- 
voltage  characteristics  were  measured  by  sweeping  the  current 
density  at  7.7  mA  cm-2  s-1  from  open-circuit  voltage  to  0.1  V  for 
each  oxygen  partial  pressure.  Resistance  was  simultaneously 
measured  at  500  Hz.  CV  was  conducted  between  0.05  and  1.2  V  at 
0.1  V  s-1  with  hydrogen  and  nitrogen  at  100%  RH  supplied  to  the 
electrodes. 

3.  Mathematical  model 

3.1.  Oxygen  transport  model 

Fig.  3  illustrates  the  oxygen  concentration  profile  across  an 
ionomer  film  on  a  planar  electrode  where  the  oxygen  reduction 
reaction  (ORR)  occurs.  When  the  diffusion-limited  current  density, 
id,  is  observed,  the  oxygen  flux  ,Jd,  can  be  described  by  Fields  law  of 
diffusion  with  constant  oxygen  diffusivity: 


Oxygen  concentration 


Fig.  2.  Schematic  illustration  of  the  cell  used  to  measure  oxygen  transport  properties  Fig.  3.  Oxygen  concentration  profile  for  an  ionomer  film  on  a  planar  Pt  electrode  when 
of  a  96  pm  thick  film.  a  diffusion-limited  current  is  observed. 
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id 

nORRp 


CSr 


Jd  =  -D02^, 


(1) 


where  Uorr  is  the  number  of  electrons  transferred  per  oxygen 
molecule  in  the  ORR ,  F  is  Faraday’s  constant,  Dq2  is  the  oxygen 
diffusion  coefficient  of  the  ionomer,  cs02  is  the  oxygen  concentration 
in  the  ionomer  at  the  gas— ionomer  interface,  and  1/  is  the  film 
thickness.  The  process  of  oxygen  dissolution  from  the  gas  phase  to 
the  ionomer  film  is  modeled  as  surface  reactions  that  are  depen¬ 
dent  on  the  concentrations  of  oxygen  in  both  phases  [30]: 


where  k02  is  a  rate  constant  and  is  the  equilibrium  concentra¬ 
tion  of  oxygen  in  the  ionomer,  which  is  a  function  of  the  partial 
pressure  of  oxygen.  From  Eqs.  (1)  and  (2), 


id  UorrFc^D q2  f  n0RRFce^2 k0l 

Therefore,  the  dissolution  rate  constant,  k0l ,  can  be  determined 
if  the  diffusion-limited  current  density  obeys  Eq.  (3)  for  various  film 
thicknesses.  An  additional  term  that  is  independent  of  film  thick¬ 
ness  appears  if  there  is  some  type  of  concentration  jump  at  the 
ionomer-electrode  interface.  The  interfacial  resistance  to  oxygen 
transport  from  ionomer-electrode  and  gas-ionomer  interfaces 
cannot  be  distinguished  from  experiments  with  a  planar  electrode. 
For  simplicity,  we  assume  that  the  resistance  exists  only  at  the  gas- 
ionomer  interface. 


3.2.  Performance  model  overview 

The  performance  model  presented  here  describes  steady-state 
species  transport  between  two  gas  diffusion  layer  (GDL)  surfaces 
of  an  MEA  at  constant  temperature.  The  MEA  is  modeled  as  a  set  of 
layers  containing  an  electrolyte  membrane,  two  catalyst  layers,  two 
microporous  layers  (MPLs),  and  two  GDLs.  The  phenomena 
described  are  (a)  the  transport  of  gas  species  (hydrogen,  oxygen, 
water  vapor,  and  nitrogen)  through  the  pores  in  the  catalyst  layers, 
MPLs,  and  GDLs,  (b)  dissolution  of  oxygen  from  the  gas  phase  to  the 
ionomer  in  the  cathode  catalyst  layer,  (c)  exchange  of  water 
between  the  ionomer  and  pores  in  the  catalyst  layers,  (d)  oxygen 
transport  through  the  ionomer  in  the  cathode  catalyst  layer,  (e) 
water  and  proton  transport  through  the  electrolyte  membrane  and 
the  ionomer  in  the  catalyst  layers,  (f)  hydrogen  oxidation  reaction 
(HOF)  in  the  anode  catalyst  layer,  and  (g)  ORR  in  the  agglomerates 
of  the  cathode  catalyst  layer.  Species  conservation  equations  in  the 
through-plane  direction  are  solved  under  given  boundary  condi¬ 
tions  at  the  GDL  surfaces  and  are  implemented  using  gPROMS 
ModelBuilder  (Process  Systems  Enterprise). 

3.3.  Conservation  law 


The  general  expression  of  the  conservation  law  for  a  steady  state 
can  be  described  as: 


dNi  * 
— -  =  0,- 
ax  ” 


(4) 


where  N*  is  the  flux  of  species  i,  x  is  the  through-plane  position,  and 
0j  is  the  volumetric  rate  of  production  of  species  i.  The  expression 
for  Nz  varies  depending  on  species  i.  The  rate  0*  is  negative  when 
species  i  is  consumed. 

The  flux  of  gas  species  through  the  pores  is  the  sum  of 
convective  flux  and  the  diffusive  flux: 


Nt  =  CXp+Ji, 


(5) 


where  C  is  the  concentration  of  the  gas,  X/  is  the  mole  fraction  of 
species  i,  v  is  the  velocity,  and  f  is  the  diffusive  flux  of  species  i.  The 
diffusive  flux  is  described  by  the  Stefan-Maxwell  equation: 


OX* 

9x 


A  XjJi-XJj 

2^  ~eff  ’ 

./  i  cd" 
j*i 


(6) 


^ejj 

where  Dy  is  the  effective  binary  diffusion  coefficient  of  species  i—j 
system  in  the  given  layer  and  M  is  the  number  of  gas  species. 

The  water  flux  through  the  ionomer  is  the  sum  of  diffusive  flux 
and  electroosmotic  flux: 


N, 


h2o 


_  _P]_neff  ^  i  4 
m7  H2°ax  +  F  : 


(7) 


where  NHl 0  is  the  water  flux,  pi  is  the  dry  density  of  the  ionomer,  mj 
is  the  equivalent  weight  (EW)  of  the  ionomer,  Dj^0  is  the  effective 
diffusion  coefficient  of  water  in  the  given  layer,  X  is  the  water 
content  of  the  ionomer,  (3  is  the  number  of  water  molecules 
accompanying  a  proton,  and  i  is  the  proton  current  density. 

The  proton  current  density  is  described  by  Ohm’s  law: 


i  =  -k 


effW 

ax’ 


(8) 


where  k^  'is  the  effective  proton  conductivity  of  the  ionomer  in  the 
given  layer,  and  is  the  electrolyte  potential.  Liu  et  al.  [31  ]  observed 
that  the  ionomer  captured  within  the  small  pores  of  primary 
carbon  particles  does  not  contribute  to  the  proton  conductivity 
across  the  electrode;  therefore,  two  models  were  compared  for  the 
effective  proton  conductivity  in  the  cathode  catalyst  layer.  One 
model,  which  will  be  referred  to  as  the  film  conduction  (FC)  model, 
assumes  that  only  the  ionomer  film  that  covers  the  agglomerate 
contributes  to  the  effective  proton  conductivity: 


£if 


1.5, 


(9) 


where  eij  is  the  volume  fraction  of  the  ionomer  film  in  the  catalyst 
layer  and  k  is  the  proton  conductivity  of  the  ionomer.  The  other 
model,  which  will  be  referred  to  as  the  bulk  conduction  (BC)  model, 
assumes  that  all  the  ionomer  in  the  catalyst  layer  contributes  to  the 
proton  conductivity: 


KeS  =  £}-5k,  (10) 

where  e\  is  the  volume  fraction  of  the  ionomer  in  the  catalyst  layer. 
Bruggeman’s  approximation  [32]  is  used  in  both  models.  The  BC 
model  is  applied  for  the  anode  catalyst  layer. 

The  HOR  rate  is  described  as  a  linearized  form  of  the  Butler— 
Volmer  equation.  The  exchange  current  density  is  assumed  to 
have  linear  dependence  on  the  hydrogen  concentration.  Lienee,  the 
volumetric  rate  of  hydrogen  production  in  the  pores  of  the  anode 
catalyst  layer  0H2,  is 


%  =  - 


s«  0  CH2  2F 
nHORFlHORcref  RuTV’ 
h2 


(11) 


where  sa  is  the  catalyst  surface  area  per  unit  volume  of  the  anode 
catalyst  layer,  nnoR  is  the  number  of  electrons  transferred  per 
hydrogen  molecule  during  the  HOR,  ij^R  is  the  reference  current 
density,  cH2  is  the  hydrogen  concentration  in  the  ionomer,  cj^  is  the 
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reference  hydrogen  concentration,  Ru  is  the  universal  gas  constant, 
T  is  the  temperature,  and  p  is  the  overpotential. 

The  volumetric  production  rate  of  water  vapor  in  the  pores  of 
the  catalyst  layer  $h2o>  is  proportional  to  the  difference  between 
the  water  activity  that  corresponds  to  the  water  content  of  the 
ionomer  a x,  and  the  water  activity  in  the  pore,  ag: 

$H20  =  Sv/<H20]^f  (aA  “  ag)>  (12) 

where  sv  is  the  area  of  the  interface  between  the  ionomer  and  the 
gas  phase,  /<h2o  is  the  evaporation  rate  constant,  and  psat  is  the 
saturation  water  vapor  pressure  at  temperature  T. 

3.4.  Agglomerate  model 

The  flooded  agglomerate  model,  in  which  the  Pt/C  agglomerate 
is  filled  with  an  ionomer,  is  applied  to  calculate  the  volumetric  rate 
of  oxygen  production  in  the  cathode  catalyst  layer,  <P02.  The  catalyst 
layer  is  comprised  of  n  spherical  particles  per  unit  volume  of  the 
catalyst  layer,  which  yields  its  porosity,  eq.  The  particle  is  comprised 
of  an  agglomerate  of  radius  R  covered  with  an  ionomer  film  of 
thickness  <5.  Hence,  n  is  given  by: 


$o2  =  n<2>!  (19) 

The  oxygen  concentration  in  the  ionomer  film  at  the  gas— 
ionomer  interface  is  treated  according  to  Eq.  (2).  The  steady-state 
flux  through  the  interface  can  be  given  as: 

Jo2  p  -l<o2(cef2  -cgf),  (20) 

where  J0l  is  the  flux  from  the  ionomer  to  the  gas  phase  and  c*+6  is 
the  oxygen  concentration  in  the  ionomer  at  the  interface.  The 
concentration  is  calculated  from  Henry’s  law.  If  the  gas— 
ionomer  interface  is  always  in  equilibrium,  then 

=  c%  (21) 

is  used  instead  of  Eq.  (20).  Eq.  (21)  is  equivalent  to  the  assumption 
made  by  Jaouen  et  al.,  that  the  oxygen  concentration  at  the  surface 
of  the  agglomerate  is  always  proportional  to  its  concentration  in 
the  gas  phase  [4]. 

For  consistency,  the  composition  of  the  cathode  catalyst  layer 
should  satisfy: 


1  -e0  =  |ir(R  +  (5)3n. 

The  area  sv  is  therefore  calculated  as: 


(13) 


sv  =  4ir(R  +  5)2n.  (14) 

Oxygen  transported  through  the  pores  in  the  catalyst  layer 
dissolves  into  the  ionomer,  diffuses  through  the  ionomer  film  and 
into  the  agglomerate  where  it  reacts.  The  species  conservation  in  an 
agglomerate  particle  covered  with  an  ionomer  film  is  described  as 
follows.  The  radial  flux  of  oxygen  in  the  ionomer  is  assumed  to  obey 
Fields  law  of  diffusion: 


j-^nar-  °^r^R 

1“°°^  R^r<R+8 


(15) 


where  J0l  is  the  oxygen  flux  in  the  radial  direction,  c0l  is  the  oxygen 
concentration  in  the  ionomer,  is  the  effective  oxygen  diffusivity 
in  the  agglomerate,  and  Dq2  is  the  oxygen  diffusivity  in  the  ionomer 
film.  The  conservation  equations  for  oxygen  are: 


neffd2C p2  _  eff  2  dCp2 

°2  df-2  02r  dr 


-  1  fef  C°2  exn 

"  ^FSaSS'0RRtf  P 

u2 

0  <  r<  R 


(16) 


n  d  Co2  n  2  dc0 

~uo2  dr2  U°?-  -  u’  - 


r  dr 


(17) 


where  sagg  is  the  catalyst  area  per  unit  volume  of  the  agglomerate, 
i'orr  is  the  reference  current  density,  cjJ  is  the  reference  oxygen 
concentration,  and  a  is  the  apparent  transfer  coefficient.  Although 
the  solution  of  Eq.  (16)  can  be  derived  analytically,  we  decided  to 
solve  this  equation  numerically  for  reasons  of  robustness. 

The  production  rate  of  oxygen  in  an  agglomerate  #1,  is  given  by: 


*1  =  4TTR2J|r_R  =  -4nR2Deff^ 


r=R 


(18) 


The  volumetric  rate  of  oxygen  production  in  the  catalyst  layer  is 
related  to  that  in  an  agglomerate  by: 


a  —  0  £o)  (  1  £fagg  y  ~Pcb  ^  £l,agg 


^ciPpt? 


(22) 


where  a  is  the  catalyst  loading,  £\<agg  is  the  volume  fraction  of  the 
ionomer  in  the  agglomerate,  y  is  I/C,  p/,  pcb,  and  ppt  are  the  true 
density  of  the  ionomer,  support  (carbon  black),  and  catalyst  (Pt), 
respectively,  and  La  is  the  thickness  of  the  catalyst  layer.  In  Eq.  (22), 
the  ionomer  volume  fraction  £iiClgg  and  ratio  <5/ft  are  left 
undetermined. 


3.5.  Parameters  for  performance  calculation 

The  following  parameters  were  used  to  calculate  the  current- 
voltage  characteristics.  The  values  given  in  Table  1  were  used  as 
the  structural  parameters  of  the  catalyst  layers.  The  ECSAs  deduced 
from  CV  were  used  for  the  cathode  catalyst  layers  (Table  1 ).  Other 
structural/physical  parameters  that  are  commonly  used  are  listed  in 
Table  2.  Thickness,  porosity,  Pt  area,  and  ionomer  parameters 
approximate  to  the  actual  parameters  for  the  MEAs.  The  rate 


Table  2 

Common  structural/physical  parameters. 


Parameter 

Value 

Thickness  of  GDL  substrate 

2  x  10"4  m 

Porosity  of  GDL  substrate 

0.8 

Anode  MPL  thickness 

6.6  x  10"5  m 

Cathode  MPL  thickness 

6  x  10-5  m 

Porosity  of  MPL 

0.5 

Membrane  thickness 

2.5  x  10-5  m 

Dry  density  of  ionomer 

1980  kg  m  3 

EW  of  ionomer 

1100  g  equiv-1 

Anode  Pt  area  per  unit  volume 

1.25  x  106  mpt  m-3 

Reference  current  density  for  HOR 

5  x  103  Ampt2 

Reference  concentration  for  HOR 

23.916  mol  m“3 

Reference  current  density  for  ORR 

5x1 0-3  A  mpt2 

Reference  concentration  for  ORR 

3.041  mol  m-3 

Apparent  transfer  coefficient  for  ORR 

1 

Rate  constant  of  water  evaporation 

100  m  s"1 
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constant  of  water  was  selected  to  be  sufficiently  large  so  that 
the  ionomer  is  almost  at  equilibrium  with  water  vapor  in  the  pore. 
The  material  properties  are  summarized  in  Appendix  A.  The 
operating  conditions  corresponded  to  those  used  for  the 
experiment. 

In  this  work,  two  models  for  the  oxygen  transport  at  the  gas- 
ionomer  interface  were  compared.  The  interface  nonequilibruim 
(IN)  model  uses  Eq.  (20),  where  the  concentration  decreases  with 
increasing  oxygen  flux  into  the  ionomer.  The  interface  equilibrium 
(IE)  model  uses  Eq.  (21)  instead  of  Eq.  (20),  where  the  oxygen 
concentration  in  the  ionomer  at  the  interface  is  always  at  the 
equilibrium  value. 


4.  Results  and  discussion 

4.1  Characterization  of  the  catalyst  layer 

Fig.  4  shows  secondary  electron  images  of  the  FIB  cross-section 
of  the  cathode  catalyst  layer  with  different  I/C,  in  addition  to  the  Pt/ 
C  (TEC10E30E)  material  (Fig.  4(b)).  Fig.  4(a)  shows  homogeneous 
structures  on  a  scale  of  ca.  1  pm.  This  homogeneity  justifies  using 
the  volume  fraction  of  the  components  that  do  not  vary  with  the 
through-plane  position  in  the  model.  The  Pt/C  particle  sizes  shown 
in  Fig.  4(b)  range  from  10  to  50  nm  or  more.  Fig.  4(c),  (e),  and  (g) 
shows  that  the  catalyst  layer  comprises  interconnecting  particles 


Fig.  4.  Secondary  electron  images  for  (a)  FIB  cross-section  of  a  cathode  catalyst  layer  with  I/C  =  0.5  at  20,000 x  magnification,  (b)  Pt/C  particles  at  100,000 x  magnification,  and 
(c-h)  FIB  cross-sections  of  the  cathode  catalyst  layer  at  100,000x  magnification.  (c,d)  I/C  =  0.5;  (e,f)  I/C  =  0.75;  (g,h)  I/C  =  1.0.  Cross-sections  are  highlighted  in  (d),  (f),  and  (h). 
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with  diameters  from  10  to  50  nm,  which  is  close  to  that  for  the  Pt/C 
particles.  These  three  images  also  show  that  voids  near  the  contact 
point  between  particles  become  filled  as  I/C  increases. 

Fig.  4(d),  (f),  and  (h)  highlights  the  cross-section  of  the  catalyst 
layer  and  reveal  that  there  are  no  large  (200-2000  nm)  ionomer- 
filled  spheres  that  correspond  to  the  agglomerate  assumed  in  the 
previous  models.  Instead,  the  Pt/C  particles  adhere  together  to 
create  the  framework  of  the  catalyst  layer  and  the  ionomer  fills  the 
voids  between  the  particles.  Although  no  large  agglomerates  were 
observed  in  the  micrographs,  relatively  tightly-packed  and  loosely- 
packed  regions  were  evident,  especially  for  large  I/C.  If  the  tightly- 
packed  region  is  regarded  as  an  agglomerate,  then  the  size  appears 
to  range  from  approximately  50-100  nm.  Although  the  size  is  close 
to  that  of  primary  Pt/C  particles,  we  assume  that  the  agglomerate 
model  is  still  applicable  to  the  experimental  catalyst  layers, 
because  a  significant  amount  of  Pt  is  supported  inside  this  type  of 
carbon  support  and  on  its  surface  [33]  and  the  Pt  utilization  is 
more  than  90%  at  80  °C  with>60%  RH  [34]. 

4.2.  Oxygen  transport  properties 

Fig.  5  shows  representative  cyclic  voltammograms  of  Pt  covered 
with  thin  Nation  films  of  different  thicknesses  under  a  diluted 
oxygen  atmosphere.  The  reduction  current  of  O2  is  observed  below 
0  V  vs.  the  pseudo-reference  electrode.  As  the  potential  decreases, 
the  reduction  current  density  initially  increases,  followed  by  an 
apparent  plateau,  and  then  increases  abruptly  with  hydrogen 
evolution.  The  reduction  current  in  the  intermediate  potential 
region  increases  with  decreasing  film  thickness.  The  concentration 


Potential  /  V  vs.  Pseudo  Ref.  Pt/02 


Potential  /  V  vs.  Pseudo  Ref.  Pt/02 


Fig.  5.  Cyclic  voltammograms  of  a  Pt  film  electrode  covered  with  thin  Nation  films  of 
various  thickness  under  (a)  1%  and  (b)  0.1%  02  at  80  °C  and  60%  RH. 


Fig.  6.  (a)  Relationship  between  the  reciprocal  of  id  and  the  ionomer  film  thickness  at 
80  °C  and  60%  RH,  and  (b)  an  enlarged  view  of  (a)  with  lines  added  to  indicate  the 
range  of  the  y-intercept. 

drop  in  the  gas  phase  is  negligible,  because  the  measured  current 
density  is  sufficiently  smaller  than  the  diffusion-controlled  current 
density  calculated  from  the  gas  diffusivity.  Although  no  clear 
current  plateau  was  observed  in  the  intermediate  potential  region, 
especially  for  thinner  films,  it  was  assumed  that  the  reaction  was 
almost  under  mass-transport  control;  therefore,  the  intersection  of 
linear  extrapolation  lines  from  the  apparent  current  plateau  and 
the  hydrogen  evolution  current  was  taken  as  the  diffusion-limited 
current,  as  indicated  by  the  dashed  lines  in  Fig.  5.  While  true 
diffusion-limited  currents  can  be  larger  than  that  estimated,  we 
assume  that  the  error  is  not  significantly  large. 


Fig.  7.  Oxygen  permeability  of  ionomer  films  at  80  °C  as  a  function  of  water  activity. 
(O):  films  with  thicknesses  between  166  and  809  nm;  (•):  96  pm  thick  film. 
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Table  3 

Oxygen  permeability  values. 


Source 

csq^Dq2  /  10"9  mol  rrr1  s"1 

T/K 

a 

Ionomer 

This  work 

0.924a 

353 

0.3 

Cast  Nafion 

This  work 

1.37a 

353 

0.6 

Cast  Nafion 

This  work 

1.96a 

353 

0.9 

Cast  Nafion 

Sethuraman  et  al.  [35] 

1.817 

353 

1 

Nafion  117 

Sethuraman  et  al.  [35] 

2.695 

348 

1 

Nafion  112 

Parthasarathy  et  al.  [36] 

0.77a 

353 

1 

Nafion  117 

a  Oxygen  solubility  is  converted  to  the  value  at  1  atm  by  dividing  the  measured 
solubility  by  the  oxygen  pressure. 


Fig.  8.  Calculated  agglomerate  radius  R,  for  the  different  models  with  I/C  =  0.5. 


The  observed  diffusion-limited  current  density  was  divided  by  the 
partial  pressure  of  oxygen  to  eliminate  the  effect  of  the  partial  pressure 
difference.  The  inverse  of  the  resulting  current  density,  i^1 ,  increases 
with  increasing  film  thickness,  as  shown  in  Fig.  6(a).  The  slope  above 
1000  nm  is  smaller  than  that  below  1000  nm.  Above  1000  nm,  hemi- 
cylindrical  diffusion  through  the  ionomer  film  can  be  considered 
comparable  to  the  linear  diffusion  normal  to  the  electrode,  of  which 
the  width  is  10  pm.  Below  1000  nm,  the  hemicylindrical  diffusion 
becomes  insignificant  and  linear  diffusion  through  the  film  is  domi¬ 
nant.  Fig.  6(b)  shows  more  detailed  results  below  1000  nm,  where 
although  the  data  are  somewhat  dispersed,  id  appears  to  obey  Eq.  (3). 

The  permeability,  c^D02,  was  obtained  from  the  slope  of  the 
regression  line  by  fitting  the  data,  and  the  results  are  shown  in 
Fig.  7  as  a  function  of  water  activity.  The  permeability  of  the  96  pm 
thick  film  is  also  shown  and  is  compared  with  the  values  in  the 
literature  in  Table  3.  The  permeability  obtained  at  a  water  activity  of 
0.9  is  between  those  reported  by  Sethuraman  et  al.  [35]  for  Nation 
112  and  Nafion  117  at  a  water  activity  of  1,  but  larger  than  that 
reported  by  Parthasarathy  et  al.  [36].  The  dependence  of  perme¬ 
ability  on  the  water  activity  is  similar  to  that  reported  by  Takamura 
et  al.  [37].  It  should  be  noted  that  the  thin  films  (<809  nm  thick) 
have  lower  permeabilities  than  thick  films.  A  similar  trend  was 
reported  by  Siroma  et  al.  [38],  where  the  proton  conductivity 
decreased  with  decreasing  film  thickness;  the  authors  pointed  out 
the  possibility  that  thin  films  absorb  less  water  than  thick  films.  If 
the  water  content  of  thin  films  is  less  than  that  of  thick  films,  it  is 
reasonable  that  oxygen  permeability,  which  decreases  with 
decreasing  water  content  [39],  is  also  depressed  in  thin  films.  We 
therefore  assumed  in  the  model  that  the  water  content  of  the 
ionomer  in  the  cathode  catalyst  layer  is  half  that  of  thick 
membranes  at  the  same  water  activity. 


S/R=  0.02 
S/R=  0.06 
S/R=  0.10 


S/R=  0.04 
5/R=  0.08 
8/R=  0.12 


Fig.  9.  Calculated  performance  curves  with  I/C  =  0.5  for  four  different  models  and  5/R  parameters.  Experimental  data  are  indicated  by  circles  for  comparison. 
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The  product  was  obtained  from  the  y-intercept  of  the 

same  regression  line  and  ranged  from  3.2  x  10-3  to 
2.0  x  10~2  mol  nrr2  s-1  at  60%  RH.  The  nonzero  value  of  the 
intercept  implies  that  the  oxygen  dissolution  process  causes 
transport  loss  in  the  cathode  catalyst  layers,  where  much  thinner 
ionomer  films  cover  the  Pt  particles.  A  comparison  of  the  first  and 
second  terms  on  the  right-hand  side  of  Eq.  (3)  shows  that  the 
dissolution  process  becomes  more  significant  than  the  diffusion 
process  when  the  film  thickness  L/,  becomes  less  than  D02/k02, 
which  is  100  nm  ( =  ce£2D02/ce£2k02  -  (lx  10-9  mol  nrr1  s_1)/(l  x 
10  2  mol  nr2  s-1))  from  the  order-of-magnitude  estimate.  The 
thickness  of  the  ionomer  covering  is  less  than  100  nm,  as  shown  in 
Fig.  4;  therefore,  the  dissolution  process  should  have  a  significant 
effect  on  the  performance  of  the  catalyst  layers.  The  resistance  to 
oxygen  dissolution,  (c^fcoj-1 ,  was  more  than  three  times  larger 
than  that  to  oxygen  diffusion,  Lf(ce£D02)-^ ,  with  a  typical  film 
thickness  in  the  catalyst  layer,  Lf  =  10  nm. 


4.3.  Agglomerate  size  and  MEA  performance 

The  agglomerate  size  cannot  be  estimated  from  the  SEM  micro¬ 
graphs  of  the  FIB  cross-section,  because  the  section  does  not  always 
pass  through  the  center  of  the  particles.  Therefore,  the  particle  size 
determined  from  the  micrographs  was  only  used  for  comparison  as 
an  alternative  method  to  determine  the  agglomerate  size.  The 
possible  range  of  8/R  can  be  determined  from  the  condition  of 
£i,agg  >  o.  For  each  8/R,  the  agglomerate  radius  R  can  be  determined  so 
that  the  model  prediction  and  the  experimental  value  are  the  same 
at  one  point  on  the  performance  curve.  The  current  density  at  ca. 
0.1  V  when  the  oxygen  partial  pressure  is  at  1.9  kPa  was  selected. 
Under  this  condition,  the  current  density  can  be  regarded  as  the 
diffusion-limited  current  density.  Fig.  8  shows  the  agglomerate 
radius  as  a  function  of  8/R  with  the  smallest  dissolution  rate, 
ceQqk02  =  3.2  x  10~3  mol  nrr2  s_1.  The  IE  model  gives  a  larger  R 
than  the  IN  model.  The  difference  between  the  proton  conduction 
models  does  not  affect  R,  except  for  low  <5/R,  where  the  FC  model 
gives  larger  ohmic  loss  than  the  BC  model.  The  value  of  8/R  was 
determined  so  that  the  model  prediction  was  close  to  the  perfor¬ 
mance  curve  when  the  cell  was  operated  at  an  oxygen  partial  pres¬ 
sure  of  19.5  kPa.  Fig.  9  shows  a  comparison  of  the  predicted 
performance  curves  for  various  8/R  with  the  experimental  data.  No 
model  prediction  fits  the  performance  curve  over  the  entire  range  of 
current  density,  but  the  models  and  8/R  parameter  can  be  selected  to 
fit  at  current  densities  less  than  0.5  A  cm-2.  The  change  in  8/R 
produces  a  small  change  in  the  performance  curve  for  the  BC  model. 
In  contrast,  the  voltage  in  the  FC  model  varies  within  the  current 
density  range,  because  the  effective  conductivity  at  the  catalyst  layer 
decreases  with  decreasing  8/R.  No  appropriate  parameter  could  be 
found  for  the  BC  model  with  I/C  =  0.5.  Therefore,  focus  was  made  on 
the  results  of  the  FC  model,  and  the  values  of  R  and  8  are  listed  in 
Table  4.  For  the  IN  model,  the  sizes  of  the  ionomer  film-covered 
agglomerates  were  from  44  to  81  nm,  which  is  similar  to  the  particle 
size  shown  in  the  micrographs  of  Fig.  4.  The  increase  in  particle  size 


Table  4 

Comparison  of  agglomerate  radius  R,  ionomer  film  thickness  5,  and  the  particle  size, 
2(R  +  8). 


MEA  Cathode  Model 


8jR  Rj nm  <5/nm  2(R  +  <5)/nm  d/R  R/nm  d/nm  2 (R  +  <5)/nm 


a 

0.5 

0.04  21 

0.9 

44 

0.08 

313 

25 

677 

b 

0.75 

0.06  27 

1.6 

56 

0.16 

210 

34 

488 

c 

1.0 

0.08  37 

3.0 

81 

0.18 

298 

54 

703 

O 

W) 

O 

> 


Current  density  /A  cm 


Fig.  10.  Comparison  of  the  predictions  obtained  using  the  IN  (solid  line)  and  IE  (dashed 
line)  models  at  an  oxygen  partial  pressure  of  9.3  kPa.  Experimental  results  are  rep¬ 
resented  by  circles. 


with  I/C  corresponds  to  the  filling  of  voids  between  Pt/C  particles.  In 
contrast,  for  the  IE  model,  the  agglomerate  particle  size  is  more  than 
400  nm  and  probably  ca.  700  nm.  However,  agglomerates  of  that 
large  size  were  not  evident  from  SEM  observation. 

Fig.  10  shows  predictions  for  current- voltage  performance  from 
the  IN  and  IE  models  with  the  experimental  data.  Although  the  IN 
model  gives  a  better  fit  to  the  experimental  results  than  the  IE 
model  for  I/C  =  0.5,  there  is  no  significant  difference  between  both 
fits.  Therefore,  both  models  can  similarly  predict  the  experimental 
results  by  the  selection  of  an  appropriate  agglomerate  size.  Thus, 
the  validity  of  the  agglomerate  model  should  be  determined  by 
how  the  model  reflects  the  actual  structure  of  the  catalyst  layer. 

From  this  discussion,  it  is  reasonable  to  attribute  the  hindrance 
of  oxygen  transport  in  the  catalyst  layer  to  slow  oxygen  dissolution 
at  the  gas-ionomer  interface,  rather  than  to  large  agglomerates. 
This  implies  that  an  increase  in  the  gas— ionomer  surface  area  and/ 
or  selection  of  an  ionomer  with  a  high  oxygen  dissolution  rate  are 
better  strategies  to  enhance  fuel  cell  performance  than  decreasing 
the  agglomerate  size. 
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5.  Conclusions 

A  model  of  oxygen  transport  in  the  catalyst  layer  of  PEFCs  was 
developed.  The  model  is  based  on  the  analysis  of  the  oxygen 
transport  properties  of  thin  ionomer  films  cast  on  a  Pt  electrode  on 
a  quartz  plate.  The  diffusion-limited  current  density  of  the  ORR  was 
measured  as  a  function  of  film  thickness.  The  current  density  could 
be  explained  by  Fields  law  of  diffusion  combined  with  the  oxygen 
dissolution  process  from  the  gas  phase  to  the  ionomer.  The  rate  of 
dissolution  was  determined  as  a  term  that  is  independent  of  film 
thickness.  The  dissolution  rate  is  significant  when  the  film  is 
thinner  than  approximately  100  nm,  from  an  order-of-magnitude 
estimate  using  the  measured  values  of  diffusivity  and  dissolution 
rate.  The  dissolution  model  was  incorporated  into  the  conventional 
agglomerate  model  with  interface  equilibrium  and  the  results  of 
these  two  models  were  compared.  The  size  of  the  ionomer  film- 
covered  agglomerate  was  determined  so  that  the  model  prediction 
of  the  diffusion-limited  current  density  corresponded  to  the 
experimental  value  and  was  compared  with  SEM  images  of  the  FIB 
cross-sections  of  catalyst  layers.  The  size  for  the  developed  IN 
model  was  close  to  that  observed  using  SEM;  however,  that  for  the 
conventional  IE  model  was  much  larger  than  the  SEM  observation. 
The  model  predictions  for  current-voltage  performance  were  not 
significantly  different.  Therefore,  it  is  reasonable  to  attribute  the 
hindrance  of  oxygen  transport  in  the  catalyst  layer  to  slow  oxygen 
dissolution  at  the  gas-ionomer  interface,  rather  than  to  large 
agglomerates. 
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(A.5) 


is  applied  instead  of  Eq.  (A.4). 

Henry’s  law  constant  for  hydrogen  in  the  ionomer,  JCH2,  is 
calculated  according  to  [43]: 

%  =  0.101325  x  0.255  x  105exp  (170/T).  (A.6) 


Henry’s  law  constant  for  oxygen  in  the  ionomer,  I<q2  ,  is  calcu¬ 
lated  as:  k0 


101325 

°2  “  4.408  -  0.09712A’ 


(A.7) 


based  on  the  data  for  a  96  pm  thick  film  cast  from  DE2020. 

The  diffusivity  of  oxygen  is  calculated  from  the  data  fit  for  the 
96  pm  thick  cast  film: 

Dq2  =  1.14698  X  10_10A°'708.  (A.8) 

The  effective  value  of  the  diffusion  coefficient  and  the  proton 
conductivity  are  described  by  Bruggeman’s  approximation  [32]. 
That  is,  the  binary  diffusion  coefficient  of  the  gas  species  in  a  porous 
component  is  given  as: 

of  =  4-5D0.  (A.9) 

The  diffusion  coefficient  of  water  in  the  ionomer  of  the  catalyst 
layer  is  given  as: 

<0  =  f5DH2o-  (A.10) 

The  diffusion  coefficient  of  oxygen  in  the  agglomerate  is  given 
as: 


Appendix  A.  Material  properties 


(A.ll) 


The  water  content  of  the  membrane  and  the  ionomer  in  the 
anode  catalyst  layer  is  calculated  according  to  the  relationship  for 
Nation  117  [38]: 

7  =  0.043  +  17.81a-  39.85a2  +  36a3  (A.l) 

However,  the  depression  of  water  content  in  the  cathode  cata¬ 
lyst  layer  is  assumed  and  approximated  by: 


7  =  1  (0.043  +  17.81a  -  39.85a2  +  36a3) .  (A.2) 

The  diffusion  coefficient  of  water  in  the  ionomer,  DH2q,  is 
calculated  from  the  self-diffusion  coefficient,  D^mr,  using  [1,40]: 


n  n  din  a 

Dr2o  =  Dnmr-^-j  exp 


2416 


J__  1 
303  ~T 


(A.3) 


The  value  of  DNMr  is  derived  from  interpolation  of  the  experi¬ 
mental  data  as  a  function  of  A  [41]. 

The  electroosmotic  coefficient  of  the  ionomer,  /?,  is  unity  for 
A  >  1.4  [42]  and  assumed  to  be  zero  at  A  =  0.  The  interpolation  is 
smooth  at  0  <  A  <  1.4. 

The  proton  conductivity  of  the  ionomer,  /c,  is  calculated 
according  to  [1] 


K 


(0.5139A  -  0.326)exp 


1268 


1 

303 


A  >  1 


(A.4) 


for  the  membrane.  The  value  at  A  =  1  is  used  for  A  <  1.  For  the 
catalyst  layers,  an  approximation  of  the  internal  data  for  a  film  cast 
from  DE2020  Nation  dispersion  solution  at  80  °C, 


Nomenclature 

a  water  activity 

C  gas  concentration,  mol  m-3 

c  concentration  in  ionomer,  mol  m-3 

D  diffusion  coefficient  in  ionomer,  m2  s-1 

D  diffusion  coefficient  of  gas  species,  m2  s-1 

F  Faraday’s  constant,  C  mol-1 

i  current  density,  A  m-2 

J  diffusive  flux,  mol  m-2  s-1 

I<  Henry’s  law  constant,  Pa  mol-1  m~3 

k  rate  constant,  m  s-1 

L  thickness,  m 

M  number  of  gas  species 

m  equivalent  weight,  kg  eq  1 

N  flux,  mol  m-2  s-1 

n  number  density  of  agglomerate,  m-3 

na  number  of  electrons  transferred  per  molecule  in  reaction 
a 

psat  saturation  water  vapor  pressure,  Pa 
R  agglomerate  radius,  m 

Ru  universal  gas  constant,  J  mol-1  K_1 

r  radial  distance  from  the  center  of  an  agglomerate,  m 

s  surface  area  per  volume,  m-1 

T  temperature,  K 

v  average  gas  velocity,  m  s_1 

X  mole  fraction 

x  through-plane  position  in  an  MEA,  m 

y  VC 


T.  Suzuki  et  al.  /  Journal  of  Power  Sources  222  (2013)  379-389 


389 


Greek 

a  apparent  transfer  coefficient 

number  of  water  molecules  accompanied  by  a  proton 
d  thickness  of  ionomer  film  over  agglomerate,  m 

e  volume  fraction 

7]  overpotential,  V 

k  proton  conductivity,  S  m_1 

A  water  content  of  ionomer 

pi  density  of  species  i,  kg  m-3 

a  catalyst  loading,  kg  m-2 

0  volumetric  rate  of  production,  mol  m3  s_1 

0  electrolyte  potential,  V 

Superscripts 

eff  effective 

eq  equilibrium 

R  at  agglomerate-ionomer  film  interface 

R  +  d  at  the  interface  between  gas  and  ionomer  film  that  covers 
an  agglomerate 
ref  reference 

s  surface 

Subscripts 
0  pore 

1  one  agglomerate 

agg  agglomerate 

CB  carbon  black 

CL  catalyst  layer 

/  film 

H2  hydrogen 

H20  water 

HOR  hydrogen  oxidation  reaction 

I  ionomer 

02  oxygen 

ORR  oxygen  reduction  reaction 

Pt  platinum 

v  ionomer-gas  interface 
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